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Abstract 



We have observed four fully reconstructed B° — > D* + D*~ candidates in 
5.8 million T(4S) — > Si? decays recorded with the CLEO detector. The 
background is estimated to be 0.31 ± 0.10 events. The probability that the 
background could produce four or more signal candidates with the observed 
distribution among D* + and D*~ decay modes is 1.1 x 10 -4 . The measured 
decay rate, B(B° -» D*+D*-) = [6.2^J(stat) ± l.O(syst)] x 10" 4 , is large 
enough for this decay mode to be of interest for the measurement of a time- 
dependent CP asymmetry. 
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The Cabibbo-suppressed decay B° — ► D* + D*~ is a promising channel for searches of 
CP violation in B° meson decays at future B factories Within the framework of the 

Standard Model, the proper time-dependent CP asymmetry in the decay B° — > D* + D*~ 
could provide a measurement of the angle (3 of the unitarity triangle |J in the same way as 
the well-known decay B° — > J/ipK° [T]J2|]. The final state D* + D*~ may be an admixture 
of CP-even and CP-odd states which could complicate such a measurement. However, the 
two CP components can be disentangled using angular correlations in the final state |4| , and 
estimates based on the heavy quark limit indicate that the dilution of the asymmetry from 
the two different CP states is small ||. The decay amplitude for the process B° — > D* + D*~ 
is expected to be dominated by the decay b — ► cW + ; W + —>■ cd. The branching fraction 
for this process can be estimated from the measured rate || of the Cabibbo-favored decay 
B° D*+D*- and is B(B° -> D*+D*-) « (f D */f D *) 2 tan 2 6 c B(B° -> D*+D*-) « 0.1%, 
where the /x are the decay constants and #c is the Cabibbo angle. 

The CLEO J7| and ALEPH @ collaborations have searched for the B° -> D* + D*- decay 
and have reported 90% CL upper limits on the branching fraction of 22 x 10~ 4 and 61 x 10~ 4 , 
respectively. In this Letter we report on the first observation of the decay B° — > D* + D*~ 
and a measurement of its decay rate. This measurement supersedes the previous CLEO 
search J7|- 

The data were recorded at the Cornell Electron Storage Ring (CESR) with two con- 
figurations of the CLEO detector, called CLEO II § and CLEO II.V. In the CLEO II.V 
configuration, the innermost wire chamber was replaced with a precision three-layer silicon 
vertex detector (SVX) [jTOf] . Each layer of the SVX is equipped with readout on both sides 
providing precise measurements of the and z coordinates of the charged particle trajectory. 
(The z-axis of the CLEO cylindrical coordinate system is coincident with the e + beam direc- 
tion.) The results in this Letter are based upon an integrated luminosity of 3.14 (2.46) fb _1 
of e + e~ data recorded at the T(4S") energy and 1.57 (1.26) fb -1 recorded 60MeV below the 
T(4S) energy with the CLEO II (CLEO II.V) configuration. The Monte Carlo simulation 



of the CLEO detector response was based upon GEANT |nj. Simulated events for the 
CLEO II and CLEO II.V configurations were processed in the same manner as the data. 

Candidates for the decay B° — ► D* + D*~ with the subsequent decays D* + — ► D°n^ and 
D* + — > -D + 7Tg were selected. The D° and D + mesons were reconstructed in the eleven decay 
modes listed in Table |. In this Letter, U D" refers to both D° and D + mesons, and "7r s " 
refers to the slow pion produced in D* + decay. In addition, reference to charge conjugate 
states is implicit unless explicitly stated. The charged track candidates from D* + and D 
meson decays were required to originate near the e + e~ interaction point. Charged kaons and 
pions were distinguished using the charged particle's measured specific ionization (dE/dx) 
and time of flight across the tracking volume. We required that the dE/dx and time-of- 
flight information was consistent with the D daughter hypotheses of the particular D meson 
decay mode. Charged tracks and K$ candidates forming a D candidate were required to 
originate from a common vertex. The K$ candidates were selected through their decay into 
7r + 7r~ mesons. The decay point of the K$ candidate was required to be displaced from the 
e + e~ interaction point and at least one daughter pion was required to be inconsistent with 
originating at the interaction point. Neutral pions were reconstructed from photon pairs 
detected in the electromagnetic calorimeter. The photons were required to have an energy of 
at least 30(50) MeV in the barrel(endcap) region, and their invariant mass was required to 
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TABLE I. The D° and D + meson decay modes used in this analysis and their branching 
fractions Q. The branching fractions B(K§ — > ir + ir~) and B{4> — > K + K~) are included for the 
modes containing K® or <j) mesons. 



D° Decay Modes 


D + Decay Modes 


Decay Mode 


Branching 
Fraction (%) 


Decay Mode 


Branching 
Fraction (%) 




3.85 ±0.09 




9.0 ±0.6 


K~TT + TT° 


13.9 ±0.9 




1.0 ±0.1 




7.6 ±0.4 


K°7r + 7r° 


3.3 ± 1.0 




1.9 ±0.1 


Kgvr+vr+vr- 


2.4 ±0.3 




3.4 ±0.4 


(j)7T + 


0.30 ± 0.03 






07T+7T 


1.1 ±0.5 


Total 


30.6 ±1.3 


Total 


17.1 ± 1.6 



be within three standard deviations of the nominal 7r° meson mass ||. The 7r° momentum 
was required to be at least 70(100) MeV for D* + (D) daughters. To reduce backgrounds, we 
accepted both (D°n+)(D°it-) and (£>%+) (L>-7T S °) combinations but not (D+ir°)(D-ir°). A 
fit constraining the mass of each D* + candidate to the nominal value improved the D* + 
momentum resolution by 14% in simulated events. 

The B° — > D* + D*~ candidates were selected by means of four observables. The first 
observable, Xmi measured the deviation of each D and D* + candidate from the nominal 
mass (M 4 n ) and mass difference (AM"), respectively, 

Xm = 

4 = 1,2 

where cx(Mj) and cr(AMj) are the average resolutions in the reconstructed D candidate mass 
Mi and the mass difference AMj = Mi(D* + ) — M iy respectively, and i = 1,2 corresponds to 
the D* + '~ and the D,D daughters. If an event had more than one B° — > D* + D*~ candidate, 
then the candidate with the lowest Xm was chosen. The second observable, L/cx(L), is the 
significance of the projected three-dimensional distance L between the reconstructed D and 
D meson decay vertices, 

\pd - m 

where pr> and Vp are the momentum and decay vertex position of the D candidate, respec- 
tively, and cr(L) was calculated from the covariance matrices of the D and D tracks resulting 
from the vertex fits of the D daughters. This observable exploits the relatively long D + meson 
lifetime and the precise decay vertex resolution available in CLEO II. V. The difference be- 
tween the energy of the B° candidate and the beam energy, AE = E(D* + ) + E(D*~) — -Ebeam, 
is the third observable. In simulated B° — > D* + D*~ decays, the AE resolution is 8.0 MeV. 
The fourth observable is the beam-constrained B° candidate mass M B = J E^ c , Am — where 
Pb is the momentum of the B candidate. The resolution of Mb, dominated by the beam 
energy spread, was measured to be 2.5 MeV [ |I~2]j . 



Mi-MfV ( AMi - AMf \ 2 
<r{Mi) ) + { a(AMi) ) 
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The selection criteria for these four observables were optimized for signal significance 
using simulated signal and background events. We also checked the optimization using 
background distributions estimated from the data by combining D candidates with 71"+ and 
7r° candidates with the momentum vector direction reversed and found similar results. The 
optimal criteria determined were x\i < 10, L/c(L) > for the (.D 7r+)(-D~7r°) candidates in 
the CLEO II. V data only, \AE\ < 20 MeV and |AM B | = \M B - M B om \ < 6.25 MeV where 
M^ om is the nominal B° meson mass [|J. 

With these criteria, the reconstruction efficiency for each D* + and D decay channel 
was measured from simulated B° — > D* + D*~ decays. Important issues in B° — > D* + D*~ 
reconstruction are the ability to reconstruct the trajectory of charged slow pions 7T+ that 
populate the momentum range from 60 to 190 MeV and the accurate determination of their 
reconstruction efficiency. The track-finding algorithm used for these results was optimized 
for the CLEO II but not the CLEO II. V configuration. As a result, the ratio of the D* + — > 
-D°7r+ reconstruction efficiency in the CLEO II. V data to that in the CLEO II data is (65 ± 
6)% due to the reduced 7r+ reconstruction efficiency. We corrected the D* + reconstruction 
efficiency for differences in the inclusive D* + meson yields between data and simulation in this 
momentum range using the measured inclusive D* + production spectrum in T(4S) — > BB 
events Including the D* + and D daughter branching fractions, the overall reconstruction 
efficiency was £ = (10.08 ± 1.10) x 10~ 4 . An appropriate figure of merit is the single 
event sensitivity, defined as [2N(BB)f 00 S]~ 1 , where N(BB) is the number of BB pairs and 
/oo = 0.48 ± O04 is the fraction of T(4S) decays to B°B° 0. Our sample of 3.3 x 10 6 
(2.5 x 10 6 ) BB pairs in the CLEO II (CLEO II. V) data gives a single event sensitivity for 
B° -> D* + D*- of (1.8 ± 0.3) x 10~ 4 . 

We used two independent methods to estimate the contributions of the background to 
the signal region, defined as |Ai?| < 20 MeV and |AM B | < 6.25 MeV. In Method 1, we 
scaled the number of candidates in a grand sideband (GSB) to estimate the background 
contribution to the signal region. The GSB is defined by the regions (50 < |AJ5| < 400 MeV 
and 5.20 < M B < 5.29 GeV) or (|A£| < 400 MeV and 5.20 < M B < 5.26 GeV) and 
is indicated in Figure [l|(a) as the area outside the dashed line. The scale factor for the 
GSB events is the ratio of the area of the signal region to area of the GSB. The estimated 
background contribution to the signal region is 0.261 ± 0.043 events from Method 1. In 
principle, the background contribution determined from the GSB slightly overestimates the 
actual background due to B — > D* + D*~X Sy d decays that are kinematically forbidden to 
populate the signal region but may be present in the AE or M B sideband regions. This 
overestimation is negligible as discussed below. 

For Method 2, we decomposed the background into four classes and estimated the con- 
tribution of each class separately. The dominant background class is composed of random 
combinations of D* + and D*~ candidates in which either one or both candidates is "fake"; 
that is, they are not composed of the daughters of an actual D* + decay. The other back- 
ground classes comprise combinations in which the D* + and D*~ candidates arise from actual 
D* + and D*~ meson decays that are roughly back-to-back. The contributing processes are 
(1) e+e" -> cc with c -> D* + and c -> D*~, (2) T(4S) -> BB with B -> D*+X and 
B — > D*~Y, and (3) B — > D* + D*~X Sy d where X s ^ represents either a strange or non-strange 
meson from the decay of an orbitally- or a radially-excited D meson or non-resonant D* + X 
production. 
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FIG. 1. (a) The AE vs the beam-constrained mass distribution for all B° — > D* + D*~ candi- 
dates in the data taken on the T(4S) resonance. The signal region is indicated by the box with the 
solid line. The area outside the dashed line is the grand sideband (GSB). There are four candidates 
in the signal region and a total of 41 candidates in the entire distribution, (b) The beam-constrained 
mass distribution for B — ► D* + D*~ candidates satisfying |A.E| < 20 MeV. 
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We estimated the combinatorial background from data with explicit fake D* + candidates 
formed by replacing Mf in \\i m equation (D) with Mf + 6a(Mi) or Mf — Qa{M\). We 
first formed a sample of fake D* + candidates combined with standard D*~ candidates. Sim- 
ilarly, we formed a sample of fake D* + candidates combined with fake D*~ candidates. The 
combinatorial background can be derived from these samples and contributes an estimated 
0.304 ± 0.040 events when scaled to the signal region. 

The contribution of the process (1) e + e~ — > cc, c — > D* + , c — > D*~ was estimated from 
the data taken 60 MeV below the T(4S) after subtracting the combinatorial background 
using the method described above. The estimated rate of e + e~ — ► cc, c — ► D* + , c — > 
was corrected for the relative cross section and luminosity and scaled to the area of the 
signal region. The estimated number of events in the signal region from this process was 
0.039 ±0.030. 

The contributions of processes (2) and (3) from T(AS) — > BB were estimated from 
a sample of simulated events approximately ten times the data sample. The process (2) 
T{AS) -> BB, B -> D* + X, B -> D*~X was estimated to contribute 0.024 ± 0.003 events 
to the signal region. The contribution of B — ► D* + D*~X S was determined to be negligible 



assuming B(B -> D* + D*~X S ) = 1.8% 0. The contribution of B -> D* + D*~X d was 



estimated from a simulation of B — > D* + D** m in which the D**^ pseudo-particle had a mass 
of 2420 MeV and a width of 400 MeV and was forced to decay to the D*~tt° final state. 
The reconstructed D* + D*~ combination from this process tends to have AE < —Mx d so the 
probability of reconstructing it with \AE\ < 400 MeV and 5.20 < M B < 5.29 GeV is twelve 
times smaller than that for the signal process. Assuming B(B — > D* + D*~Xd) pa B(B° — > 
D* + D*~), this contribution to the signal region is negligible. The estimated contribution to 
the signal region from the sum of all backgrounds is 0.367 ± 0.051 events for Method 2. 

The background rates obtained from these two statistically independent methods were 
averaged to yield the estimated background contribution to the signal region of 0.306 ± 
0.033(stat) ± 0.094(syst). The 31% systematic uncertainty arises from the uncertainty in 
the shapes of the AE and Mb distributions of the background. The systematic uncertainty 
was taken to be the difference in the scale factor when these distributions were fitted with 
second- and first-order polynomials, respectively, instead of a zeroth-order polynomial. 

The distribution of the 41 candidates passing the selection criteria in the T(45) data 
sample in the AE vs Mb plane is shown in Figure [TJ. There are four candidates in the 
signal region. The observed number of candidates and the estimated background for the 
(L) 7r s + )(U 7r & 7) and (D°7i+)(D-^) submodes are listed in Table |TJ. Also listed in Table is 
the probability that a fluctuation of the estimated background could produce the observed 
number of signal candidates or more in each submode. The calculation of the background 
fluctuation probability assumes that the statistical uncertainty in the background in the two 
submodes is uncorrelated and that the systematic uncertainty in the background is fully 
correlated between the submodes. Integrating over all background levels, assuming that 
the number of background events is normally distributed about its central value for each 
submode [IB], we find that the combined probability that the estimated background could 
produce the observed number of signal candidates or more in the two submodes is 1.1 x 10~ 4 . 

The branching fraction of -B — > D* + D*~ was calculated using a maximum likelihood 
technique that took into account the signal efficiency and estimated background contribution 
for each D decay mode in the CLEO II and CLEO II. V data samples. The branching fraction 
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TABLE II. The efficiency, observed number of candidates and estimated number of background 
events in the (D°ir+)(D°ir~) and (L>°7r+)(-D _ 7Tg ) decay submodes. The reconstruction efficiency 
£ includes the D and D* + daughter branching fractions. The row labeled "All (AE,Mb)" is the 
total number of B° -> D*+D*~ candidates in each submode in the 5.20 < M(B) < 5.29 GeV 
and \AE\ < 400 MeV region. The row labeled "Signal region" contains the observed number of 
signal candidates in the \AE\ < 20 MeV and \AM B \ < 6.25 MeV region. "Bkg. Meth. 1" and 
"Bkg. Meth. 2" are the number of background events in the signal region estimated using the two 
independent methods described in the text. The sixth row contains the average estimated number 
of background events in the signal region. Only statistical uncertainties are included for the upper 
six rows. The calculation of the background fluctuation probability P is described in the text. 









Total 


£ x 10 4 


6.06 ± 1.02 


4.02 ± 0.40 


10.08 ± 1.10 


A\l(AE,M B ) 
Signal Region 


13 
2 


28 
2 


41 

4 


Bkg. Meth. 1 
Bkg. Meth. 2 


0.080 ± 0.024 
0.091 ±0.024 


0.181 ±0.036 
0.275 ± 0.044 


0.261 ± 0.043 
0.367 ±0.051 


Average Bkg. 


0.085 ± 0.017 


0.219 ± 0.028 


0.306 ± 0.033 


P 


3.85 x 10" 3 


2.24 x 10~ 2 


1.10 x 10~ 4 



was determined to be B(B° -> D*+D*~) = (6.2±|;g ± 1.0) x 10" 4 , where the first error is 
statistical and the second is systematic. The systematic uncertainty is dominated by the 
uncertainties in S (13.3%) and /oo (8.3%). The product branching fraction, using the D* + 
and D decay modes as in this Letter, B(B° -> D* + D*-) x ££(£>*+ -> Dir) x B(D -> 
X) 4 x 10" 5 , is comparable to that of the B° — > J/ipK° decay mode, B(B° — > J/ipK°) x 
EB(J/i/j -> £+£-) x B(K° -> ^ -> 7T+7T-) w 3.6 x 10~ 5 . With careful optimization of the 
charged track reconstruction efficiency, in particular that of charged slow pions, exclusively 
reconstructed B° — > D* + D*~ decays could permit a complementary measurement of the 
angle (3 of the unitarity triangle at future B factories. 

In conclusion, we have fully reconstructed four B° — > D* + D*~ candidates with a total 
estimated background of 0.31±0.10 events in 5.8 x 10 6 T(4S*) — > BB decays. The probability 
that the estimated background could fluctuate to the observed number of signal candidate 
events or more is 1.1 x 10~ 4 . The branching fraction is measured to be B(B° — > D* + D*~) = 
[6.2^(stat) ± l.O(syst)] x 10~ 4 and is in agreement with the expected rate. This rate 
suggests that this decay could provide an avenue for the measurement of the angle j3 of the 
unitarity triangle. 
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